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Three-body decay functions in time-like parton branching are calculated using the 
jet calculus to the next-to-leading logarithmic (NLL) order in perturbative quantum 
chromodynamics (QCD). The phase space contributions from each of the ladder dia- 
grams and interference diagrams are presented. We correct part of the results for the 
three-body decay functions calculated previously by two groups. Employing our new 
results, the properties of the three-body decay functions in the regions of soft partons 
are examined numerically. Furthermore, we examine the contribution of the three-body 
decay functions modified by the restriction resulting from the kinematical boundary of 
the phase space for two-body decay in the parton shower model. This restriction leads 
to some problems for the parton shower model. For this reason, we propose a new 
restriction introduced by the kinematical boundary of the phase space for two-body 
decay. 
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1. Introduction 



The next generation of linear colliders (NLC), such as the Japan linear collider (JLC) ^ 
and the super conducting electron positron linear collider (TESLA) 2 ) will be running in the 
near future. The energies of these accelerators will be in the TeV region. Many multijets 
will be produced in this energy region. In order to analyze the experimental data for such 
multijet events, parton shower models based on perturbative QCD will be useful. To this 
time, many parton shower models limited to the contribution of the leading logarithmic (LL) 
order have been constructed. 3 ) 

With the purpose of constructing a parton shower model including the contribution 
of the next-to-leading logarithmic order, three-body decay functions have been calculated 
previously by two groups. 4 )' 5 ) Using these results, NLLJET has been constructed, and it 
has been used to analyze experimental data, such as that from SLAC-PEP, TRISTAN and 
LEP1. 6 ) 

At present, only two groups have calculated the three-body decay functions. Since some 
discrepancies have been found in the results reported in Refs. 4) and 5), the results for the 
three-body decay functions must be verified. 

The three-body decay functions in the parton shower model have been modified with 
respect to the manner in which the effect of the one soft gluon is absorbed into the kinematical 
boundary in the phase space of the two-body decay. 9 ^ This modification corresponds to 
the angular ordering of the three-body decay functions. 

We expect that many soft partons will be produced in the TeV energy region. In order 
to decrease the systematic error involved in the data analyses of multijet events, the effects 
of the soft partons should be included in the parton shower model. For example, there is 
the double soft gluon limit. Theoretical studies of the soft parton for the three-body decay 
functions are necessary to improve the parton shower model, such as NLLJET. 

In this paper, we employ the method of Ref. 7) to construct the kinematics of three-body 
decay functions in time-like branching. The contributions of the three-body decay functions 
for the soft partons are then examined using the results so obtained. 

The organization of the paper is as follows. In the next section, we construct the phase 
volume element. The details of the kinematics are given in Appendix A. In §3, the calculated 
results are presented. We find that part of the results obtained in this paper do not agree 
with those obtained previously by two groups. In §4, the properties of the calculated three- 
body decay functions are examined. In §5, we give some comments on the angular ordering. 
In the final section, we give a summary of our conclusions. 



2 



§2. Kinematics for the calculation 

In this section, the kinematics used to calculate the three-body decay functions are 
presented. The three-body decay functions are determined from the time-like jet calculus 
on the basis of the QCD factorization theorem. 4 )' 5 )' 10 ) To begin with, the factorization is 
introduced following Ref. 10). In the light-cone gauge, collinear singularities appear from 
the on-shell partons (parton internal lines) between two-particle irreducible (2PI) kernels, 
and they are factorized using the projection operator. 

We use A and B to represent the 2PI kernels that are connected by gluon internal lines 
or quark internal lines. The amplitude squared C is defined as 

C = VS" V (2-1) 

for gluon internal lines and 



C = Tr[AB] 



(2-2) 



for quark internal lines, where // and v' are Lorentz indices. The projection operators to 
extract the collinear singularities are obtained from the general formulas of the 2PI kernels. 10 ^ 
We use q to represent momenta on the internal lines that correspond to the on-shell partons. 
Inserting the projection operator p( G \q) for gluon internal lines into Eq. (2-1), we have 



P (G \q)C 



A^d^'iq) 



q 2 =0 



(2-3) 



where \i and v are Lorentz indices, the factor 1/2 is the average over gluon helicities, and 
d^' v ' ' (q) is given by 



d^(q) = -g^ + (r// + q»'n v ')l{nq). 

Here n is a light-like vector that specifies the light-cone gauge. Similarly, inserting the 
projection operator P^\q) for quark internal lines into Eq. (2-2), we have 



P^\q)C = Tr 



A aa '{ 



__ Q Aqn 



-Tr 



(2-4) 



where a, (3, a' and f3' are the spinor indices. Hence, using Eqs. (2-3) and (2-4), the projection 
operators are given by 
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(2-5) 



and 



(2-6) 



The second factors on the right-hand sides of Eqs. (2-3) and (2-4) include the collinear 
singularities. 

The diagrams contributing to the corrections at order a 2 s in the light-cone gauge are 
displayed in Fig. 1. 

In the jet calculus, 5 - ) the minimum mass scale Mq is used as the resolution that distin- 
guishes one jet from two jets in the connected line for the ladder diagram. For instance, Mq 
can be applied to the lines labeled "j" and "k" in Fig. 1. The three-body decay function 
arises from the contribution above Mq. Also, the contribution below Mq is absorbed into 
the two-body decay. 5 ) 




+cyclic 



(C) 












k 


+cyclic 



Fig. 1. The diagrams contributing to the corrections at a? order in the light-cone gauge. 
The labels i, j and k represent the momenta of the partons. (A) Type [L]: Ladder di- 
agrams with two identical time-like propagators. Here, L^cxl/t^, L^ocl/tfg, L^ot\/t\ 2 - 
(B) Type [I]: Interference diagrams with two different time-like propagators. Here, 
I^o<:l/ti2t23, /( 2 )ocl/ii2ii3, 1^ oc l/ti3t23- (C) Four-gluon interaction for pure gluon decays. 

In this paper, we calculate the three-body decay functions for the following processes: 
Process(Gl) 

g(q)^g{h) + g(k) + g(h), 
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Process(Ql) 

q(q)^g(li)+g(l 2 )+q(k), 

Process(G2) 

g(q)^q(h) + g(k) + q(k), 

Process(Q2) 

q(q)^q(h) + q(k) + q(h)- 

Here, the momentum of the mother parton with time-like virtuality is denoted by q. The 
momenta of the three daughter partons on the mass shell are denoted by l±, I2 and Z 3 . The 
kinematics for the variables are determined as q 2 = t and — — l\=^- 
Here, the invariants are defined as 



Uj = (li + lj) 2 , (¥J) (2-7) 
ti2 + t n + t 2 s = t. (2-8) 

The collinear contribution of the branching vertex at a 2 order extracted by the projection 
operator is defined as 



V = g A JdrP^M^ (2-9) 



where g is the QCD coupling constant, M is the squared matrix element, and [where 
i=Q (quark) or G (gluon)] is the projection operator that is given by Eqs. (2-5) and (2-6). 
Also here, dr is the three-body phase space volume element, which is given by 



The Sudakov variable for the final state parton is defined as 



U = ann + Pin + l iT , (n 2 = n 2 = nliT = nl iT = 0) (2-11) 
where n is a vector, Ut is the transverse momentum, and 



a,- 



J 2 

1 iT 

2f3ifin 



(2-12) 
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The Sudakov variable for the parton with virtuality t is also given by 



q = a q n + f3 q n + q T , (2- 13) 

where qT is the transverse momentum and we set qT — 0, a q — f3 q — 1. 
Using Eqs. (2-11) and (2-13), the momentum fraction Zi is defined by 



Zi = — = — = A- (2-14) 
qn nn 

From the Sudakov variables and the momentum fractions, the phase space volume element 
can be written 



dr = ^-y 6 5(q - k - l 3 - l k )^^d% T d%Td%5(ll), (2-15) 

with k) = (1, 2, 3) or (1,3,2). This phase space volume element can be written in terms 
of dr, which is defined as 



df = d 2 hTd 2 l jT S(ll), (2-16) 
with (i,j,k) = (1,2,3) or (1,3,2). Employing Eq. (2-16), we find that dr takes the form 



dr = 1 - Z-L-Z2- z 3 )dzidz 2 dz 3 dr, (2-17) 

where X = ZiZj and = (1,2) or (1,3). 

From Eq. (2-17), the branching vertex function is given by 



(\ 2 , , 

^ J 5{l-z 1 -z 2 - z 3 )d Zl dz 2 dz 3 Jj, (2-18) 

where a s = (? 2 /47r, and 



J =^x dtp&M h- (249) 

The quantity J is called the "decay function." It consists of a single logarithmic term and 
a double logarithmic term, due to the collinear singularity. The contribution of the double 
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logarithmic term should be subtracted from that of the three-body decay function to avoid 
double counting. 4 )' 5 )' 10 ) In order to determine the distribution for 5(1 — z\ — z<i — zs) dz\dzi dz% 
in Eq. (2-18), we integrate over the phase space with volume element dr. 

We now briefly explain the calculation of the squared matrix element P^M. We used 
REDUCE *) to calculate this squared matrix element. Replacing the inner product of the 
momenta in P®M by the momentum fractions and the invariants, P^M can be calculated. 
[See Eqs. (A-6) and (A-23) in Appendix A.] 

Integrating P^M over the azimuthal angle, it can be verified that there are no terms 
proportional to t 2 in the numerator of P^M. [See Eq. (A- 11) in Appendix A.] The ladder 
diagrams cannot have terms that are more singular than logarithmic order. Furthermore, the 
condition that the interference diagram has no mass singularity was confirmed. We examined 
whether or not the terms proportional to the invariants in the numerator of P^M satisfy 
the condition of having no mass singularity. [See Eq. (A-32) in Appendix A.] The details of 
the calculations are presented in Appendix A. 



§3. Calculated results 

Here the calculated results for the decay functions in the light-cone gauge are presented. 
The ladder diagrams include the terms which are the convolutions of the Altarelli-Parisi 
splitting functions. u ) Each type of calculated decay function is represented as 



jm= f L f^ + L f> l0gVt + Lf, (3-1) 

J M 2 tij 

j[^W?W^ + 4 fe) , (i,j,k = l,2,3,i^k) (3-2) 

Zk 

where yi — 1 — Zi (i = 1,2,3) and 0(MQ/t) terms are ignored. Mq is the minimum mass 
scale of the phase space integration. 

Integrating over for Type [L] and summing over all contributions from J^ L(1 ^ to j' /(3) l, 
we obtain 7 ^ 

/(3) 

£ J^ = V LL \og-t + V NLL , (3-3) 

i=LW M 

where 

Vll = L ( P + 4 2) + Lf (3-4) 

*) REDUCE was developed by A. Hearn. 
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and 

Vnll = V l + V n , (3-5) 

with 

V L = j:Lf]D S y i + £ J^log^, (3-6) 

i=l j,k,l=l z i 

^ = E4' + E#. (3-7) 

i=l j=l 

The terms Vll and Vnll m Eq. (3-3) are the contributions to LL order and NLL order at 
order a 2 sl respectively. The term Vnll-, called the three-body decay function, is constructed 
from the logarithmic term Vl and the non- logarithmic term Vn- 

The three-body decay functions are calculated in Refs. 4) and 5). We compare the results 
given there with the new results reported in this paper. The results for Processes (Gl) and 
(Q2) in this paper are consistent with the calculated results in Refs. 4) and 5). While the 
results for Process (Ql) obtained in Refs. 4) and 5) are quite different, our result agrees with 
that in Ref. 4). The results, excluding the logarithmic terms for J*- 1 -* and I^, for Process 
(G2) obtained in this paper are consistent with those in both Refs. 4) and 5). 

Using the notation defined in §2 and 3, the results for the logarithmic terms of the 
interference diagrams 1^ and for Process (G2) are presented. The results used in 
NLL JET 4 ) are 



/W log™ +/w i og 



2/22/3 



(3) 



ym 



—CaTr 
— CaTr 



4 + 4 

z 2 



z 5 

v2 i ^2 



zf + zj + 1 

2/2 

z\ + z\ + 1 



+ 1 - Zi + z 3 



log 



2/22/3 



Z2 2/2 

For comparison, the results given in Ref. 5) are 



+ 1+^1-^3 



log 



Z3 

ym 

Zl 



I? ^g 



Zl 



zz 

2 i Jl 



— —CaTr 
— CaTr 



zf + z^ zf + zi + 1 



z 2 



2/2 

v2 I ^2 



+ 1 + Zi - Z 3 



log 



2/22/3 



Zl 



zj + zj + 1 

2/2 



+ 1 - Zi + z 3 



log 



Zl 

2/12/2 

Z3 



(3-8) 



(3-9) 
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The new results obtained in this paper are the following: *) , 8 ) 



2/22/3 



r (i) , </2'</3 . r (3) , 
I) log — — h r L log 



2/12/2 



—CaTr 
— CaTr 



Zl 
' 2 

Z l 



Z 3 



J2 
z 3 



2 

Z l 



z 3 



1 



Z2 

z\ + zl 



2/2 

z\ + z\ + 1 



1 - 



^2 



2/2 



+ 1 + 



Zl - z 3 

2/2 

2i - z 3 



log 



2/22/3 



2/2 



log 



Zl 

ym 

Z3 



(3-10) 



0.3 



0.2 



10" 



-10" 



-0.2 
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Fig. 2. The Z2 dependence of R described by Eq. (3T1) at z% = 10 1 for Process (G2). The 
symbols correspond to the different numerators for R corresponding to Eqs. (3-8)~(3-10). 

We next examine the numerical results for the calculated three-body decay functions 
given in Refs. 4) and 5) and the new results of this paper. We define the ratio R as 



jW log ym + jw log m 



r(3) 



R 



LL 



(3-11) 



Equations (3-8)-(3-10) are inserted into the numerator in Eq. (3-11). As shown in Fig. 2, 
the numerical result for Eq. (3-10) differs significantly from the results given in Refs. 4) and 
5) in the region [-2 2 >1CT 2 ] for z 3 = 1CT 1 . 



§4. Numerical results 

In order to investigate the properties of the three-body decay functions in the soft re- 
gions, some numerical quantities are examined. These properties are useful to construct and 

*) As discussed in Ref. 7), using the results given in Ref. 8), only Eq. (3-10) satisfies the relation of the 
crossing symmetry between the time-like branching process and the space-like branching process. 
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improve the parton shower model. 
We define the following ratios: 



Ri = 




(4-1) 



R2 



V L 



(4-2) 



(4-3) 



If the absolute value of R\ is large, the contribution of the NLL order becomes important, 
because the contribution of the term Vnll is large relative to that of the term Vll at a 1 
order. If Ri~R 2 , the dominant contribution of Vnll is the logarithmic term Vl- The ratio 
of the contribution of the logarithmic term Vl to that of the non- logarithmic term Vn is 
represented by -R3. 

4.1. z 2 dependences of Ri and R 2 

In this subsection, we focus on the relation between the contributions of the logarithmic 
term Vl and the non- logarithmic term Vn- In Figs. 3-5, the gluon momentum fraction z 2 
dependences of R\ and R 2 are presented at z 3 = 0.5, 10 _1 , 10~ 2 , 10~ 3 . Here, R\ and R 2 are 
depicted by the solid curves and the crosses, respectively. 

First, we examine the z 2 dependences of R\ and R 2 for Process (Gl), in which a gluon 
decays into three gluons. In Fig. 3, R\~R 2 in most of region depicted. The absolute values 
of Ri and R 2 in the double soft gluon region [z 2 — ^3— yi/2<0] are small relative to those in 
the other region. 

The numerical result for Process (Ql), in which a quark decays into one quark and two 
gluons, is shown in Fig. 4. Here, R\~R 2 , except in the region of the small quark momentum 
fraction Z3 and the hard gluon momentum fraction z 2 (soft quark region). The contribution of 
the non- logarithmic term Vn is important in the soft quark region. As the quark momentum 
fraction z 3 is hard and the gluon momentum fraction z 2 is small (soft gluon region), the 
absolute values of R\ and R 2 are large. In addition, the absolute values of R\ and R 2 in the 
soft quark region are small. 

The numerical result for Process (G2), in which a gluon decays into one quark, one 
anti-quark and one gluon, is presented in Fig. 5. In Process (G2), Ri~R 2 in the strong 
ordering region [z 2 <^zs<^zi\. As the anti-quark momentum fraction z 3 is soft and the gluon 
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-3.5 



10"' 



10 ; 



icr : 



10- 



Fig. 3. The z<i dependences of the ratios R\ and R 2 for Process (Gl). The solid curves and the 
crosses denote R\ and R 2 - The symbols □, *, x and + correspond to the values 1CT 3 , 10~ 2 , 
10 _1 , 0.5 for z 3 . 



momentum fraction z 2 is hard (soft anti-quark region), the absolute value of Ri is large 
relative to that of R 2 . Thus, the contribution of the non- logarithmic term gives an important 
contribution to the soft anti-quark region. 

In the following subsections, we examine the properties of the three-body decay functions 
in each of the R\~R 2 regions and the soft fermion (quark or anti-quark) regions. 




11 



0.5 


-0.5 
-1 

"i -1.5 

cc 

-2 
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-3.5 

10' 4 10~ 3 10" 2 10" 1 1 

Z 2 

Fig. 5. The Z2 dependences of the ratios R\ and R2 for Process (G2). The solid curves and the 
crosses denote R\ and R2. The symbols □, *, x and + correspond to the values 10~ 3 , 10~ 2 , 
10 _1 , 0.5 for z 3 . 

4.2. Cancellation of the non-logarithmic term 

As in Ref. 7), we examine the cancellations of the non-logarithmic terms VJv in detail. 
Figures 6-8 confirm the structure of the cancellations of the non- logarithmic terms Vn in 
the region in which i?i~i? 2 . For each process, the structure of the cancellation of the non- 
logarithmic term V N is different. For these processes we have the following relations: 

Process (Gl) 

(7« + 4 2 ) + 4 3 XLff + Lg)+Lg)), (44) 

Process (Ql) 

(lV + lSy(LW + L% + Lf+lP), (4-5) 

Process (G2) 

(/<; , + /}JXLS;» + Li?' + Lg , + «»). (4-6) 

In the case of Process (Gl), Eq. (4-4) implies that the contributions for the non- 
logarithmic terms of the interference diagrams are canceled by those of the ladder diagrams. 
Equations (4-5) and (4-6) imply that the contributions of the three ladder diagrams and one 
interference diagram for Processes (Ql) and (G2) cancel those of the non-logarithmic terms 
for the interference diagrams. 
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Fig. 6. The Z2 dependences of non-logarithmic terms for Process (Gl) at zs = 10 1 . The notation 
(Ln, etc.) is defined by Eq. (3-7) in the text. Here, Z2=10 _4 -0.9. 
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Fig. 7. The Z2 dependences of non-logarithmic terms for Process (Ql) at z% = 10 1 . The notation 
(Ljv, etc.) is defined by Eq. (3-7) in the text. Here, Z2=10 _4 -0.9. 



4.3. Soft fermion [quark-anti- quark] region 

As mentioned above, we showed that the dominant contribution of Vmll is the logarithmic 
term Vl, except in the soft quark region for Process (Ql) and the soft anti-quark region for 
Process (G2). First, we examine the soft quark limit for Process (Ql). In Fig. 9, there are 
some peaks of the ratio R 3 present in the soft quark region [-2 3 <1CT 4 ]. In order to understand 
the meaning of the peaks in Fig. 9, we examine the structures of the logarithmic terms Vl 
in the soft quark region in detail. 

In Fig. 10, the contribution of the logarithmic term V L for the soft quark region crosses 
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Fig. 8. The z 2 dependences of non-logarithmic terms for Process (G2) at z 3 = 10" 1 . The notation 
(Ljv, etc.) is defined by Eq. (3-7) in the text. Here,Z2=10 -4 -0.9. 
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Fig. 9. The z 2 dependence of the ratio R3 with z% = 10 4 , 10 5 and 10 6 for Process (Ql). Here, 
z 2 =0.05-0.8. 

zero at some points. As the quark momentum fraction z 3 is small, the contribution of the 
logarithmic term crosses zero in the small z 2 region. 

Next, we show how the logarithmic term Vj, crosses zero. The momentum fractions in 
the soft quark region are given by 



J/3-1, 2/1—^2, V2^Z 1 . 

From Eq. (4-7), the logarithmic term V L can be approximated as 



(4-7) 



14 




Fig. 10. The z 2 dependence of V L for Process (Ql) with z 3 = 10~ 4 , 10~ 5 and 1(T 6 . Here, z 2 =0.05- 
0.8. 



with 



VL-L^tog^ + Lflogz, + /flog^, (4-8) 

^3 



L^Cl^L^C* 1 -^ (4-9) 



and 



j£ 3) ~ - 2CV(c F - l -C A )—. (4-10) 
Using Eqs. (4-9) and (4-10), the dominant singular term in Eq. (4-8) can be given as 

Vl~ log — — H log . (4-11) 

z 2 y 2 ziz 2 z 2 y 2 z 3 

In Fig. 11, since the contributions for these two terms cancel, the logarithmic term Vj, 
crosses zero. The first term and the second term in Eq. (4-11) arise from the ladder diagrams 
and the interference diagram, respectively. 

As seen from Fig. 12, the peaks of the ratios R 3 for Process (G2) appear near z 2 = 0.6. 
As for Process (Ql), it is expected that the contributions of the logarithmic terms are small 
in this region. 
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Fig. 11. The z 2 dependences of the logarithmic terms of the ladder diagrams and the inter- 
ference diagrams /(*) for Process (Ql) with Z3 = 10 5 . The logarithmic terms for the ladder 
diagrams L^ k \ etc. are given by Eq. (3-6). Here, Z2=10 _1 -0.8. 




10- 1 

z 2 

Fig. 12. The z 2 dependence of the ratio R 3 for Process (G2) with z 3 = 1(T 3 , 1(T 4 , 10" 5 and 1(T 6 . 
Here, z 2 =0.05-0.8. 

In Fig. 13, because the sign of the numerical value of the logarithmic term for the ladder 
diagram of Process (G2) is opposite to that of the sum of the logarithmic terms for the 
interference diagrams 1^ and I^, the contribution of the logarithmic term Vl is small. It 
is easy to derive this structure analytically. The logarithmic term Vj, in Eq. (4-8), taking 
account of the limit z 2 ^l and the relation z^<^z 2) is given by 

V L ~L { *hogy 2 + I<j?hog^, (4-12) 

Z3 
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with 



vl 



2z 2 -l 

Z2V2 



(4-13) 



In Eq. (4-13), the contribution of the logarithmic term for the interference diagram J® is 
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Fig. 13. The z 2 dependences of the logarithmic terms for the ladder diagrams and the inter- 
ference diagrams /(*) for Process (G2) with z% = 10 5 . The logarithmic terms for L^ k \ etc, are 
given by Eq. (3-6). Here, z 2 =10 _1 -0.8. 

positive for z 2 >0.5, which agrees with the behavior depicted in Fig. 13. 

In the next section, using the results obtained in this section, we discuss how to treat 
the contributions of the three-body decay functions in the parton shower model. 



§5. Comments on the angular ordering 

5.1. Angular ordering for the soft gluon region 

Before the modifications of the three-body decay functions in the parton shower model 
are discussed, we fit the approximated forms to the numerical results in order to check the 
results obtained in §4. In the region of the soft gluon momentum fraction z 2 , the momentum 
fractions are given by 

yi^z 3 ,y 2 ^l,y 3 ~z 1 . (5-1) 
By use of Eq. (5-1), the logarithmic term Vl is approximated by 



17 



with 



VL-L^logzn + Lf\og Zl + if log 



(5-2) 




(5-3) 



for Process (Gl), 




(54) 



for Process (Ql), and 




f^2C F T R H( Zl )-, -lf~4T R (C F - l -C A )H( Zl )- 



(5-5) 



for Process (G2), where P{z) = (1 + ^ 2 )/(l - z),K(z) = 1/z + 1/(1 - z) - 2 + z{\ - z) 
and H(z) = z 2 + (1 — z) 2 are Altarelli-Parisi splitting functions. 11 -* The contributions to 
Eqs. (5-3)-(5-5) agree with the solid curves R\ in Figs. 2-4. Also, since the contributions of 
the non-logarithmic terms Vjv in the one soft gluon regions are small relative to those of the 
logarithmic terms V F (Figs. 2-4), we ignore the contributions of the non- logarithmic terms 



Equations (5T)-(5-4) offer a physical interpretation of the contributions of the three-body 
decay functions for the parton shower model. The singular factor log zjz of the three-body 
decay function in the one soft gluon limit comes from the contribution of the interference di- 
agram. This contribution due to the interference diagram causes the perturbative expansion 
to break down. In addition, as shown in Figs. 3 and 4, the contributions of the three-body 
decay functions are negative for small z 2 for Processes (Gl) and (Ql). Since the parton 
shower model generates momentum fractions of the branching partons by decay functions, 
the contributions of the decay functions are required to be positive. In order to avoid these 
difficulties, it is necessary for the singular term to be absorbed into the kinematical boundary 
of the phase space for the two-body decay by the angular ordering. 4 )' 9 ) 

We expect that many soft gluons will be produced in the TeV energy region. There 
should exist the situation that the momentum fractions of the two gluons are soft. For 
this reason, we examine the relation between the double soft gluon limit and the angular 



V N . 
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ordering reported in Ref. 4). We obtain momentum fractions for the double soft gluon limit 
that satisfy the following: 



^i~z 2 ~y, z 3 ~l. (5-6) 

As shown in Figs. 3 and 4, R 1 ^R 2 in the double soft gluon region. The dominant contribu- 
tion of Vnll m this region is the logarithmic term V L . It yields 



\4~4 3) logy 3 . (5-7) 

Therefore, we calculate the term in the double soft gluon limit for Process (Gl) and 
(Ql). We find 



2/3 

for Process(Gl) and 



T^C\±- 2 (5-8) 



L fJ-C F C A P(z,)- (5-9) 
2 2/3 

for Process(Ql). The angular ordering in the double soft limit reported in Ref. 4) satisfies 
these approximate formulas for the double soft gluon limit, Eqs. (5-8) and (5-9). 

5.2. New restriction introduced by the kinematical boundary in the phase space for Process 
(G2) 

As shown in Fig. 5, the contribution of the three-body decay function for Process (G2) 
is negative in the soft anti-quark region. For Process (G2), both the logarithmic terms 
and the non- logarithmic terms for the ladder diagrams were absorbed into the kinematical 
boundaries of the phase spaces for the two-body decays to maintain the positivity of the 
three-body decay function. 4 ) 

Following Ref. 4), the formulas of the non-logarithmic terms for the ladder diagrams are 
written 



-(*) 



-4 fc) + u. 
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(5-10) 



Therefore, a part of the non-logarithmic terms can be absorbed into the kinematical bound- 
aries of the phase spaces for the two-body decays. By use of Eq. (5-10), we obtain 



f Lni _^ + L , hogVt + L , 



" "*J 

ykt L^(l-^)^l + L>, (5-11) 

with = (1,2), (1,3) or (2,3). The contribution of 0{Mq/\) is absorbed into the two- 
body decay. 

We corrected a part of the results for Process (G2) in §3. The logarithmic terms for the 
interference diagrams 1^ and J*- 3 -* used in NLLJET are replaced by Eq. (3-10). Now, we 
need to examine the contribution of the three-body decay function for Process (G2) given 
in Ref. 4) modified by Eq. (5-11). It is given by 

V NLL - J2[ L l } l °SVi + L { $(CV) 
W m = ^ , (5-12a) 



with 



V, 



£l$(CV) = -C F T R H( yi )P(^)-- C F T R H(y 3 )p(^)- 

f^i K yi J yi K ys ; y3 

-2C A T R K(z 2 )H(^)-, (5-12b) 
v 2/2 y 2/2 



where L^(CV) are the convolutions of the Altarelli-Parisi splitting functions u ) in the non- 
logarithmic terms for the ladder diagrams. We call the function in Eq. (5- 12a) the modified 
three-body decay function. We compare this decay function with the LL term in order to 
explore the contribution at the three-body level. Both the first term and the second term in 
Eq. (5- 12b) differ from the corresponding results in Ref. 4). The results in Ref. 4) are not 
the convolutions of the Altarelli-Parisi splitting functions. They are given by 



CfTr 



(5-13) 



ff W p(i)! + fl(„)P(*i)± 

Although the numerical difference between Eq. (5-13) and the convolution term in Eq. 
(5-12b) is small, we adopt the result Eq. (5-12b). 

As shown in Fig. 14, substituting Eq. (3-10) into Eq. (5-12a), the contribution of Eq. 
(5T2a) is negative for the anti-quark momentum fraction z 3 = 1CT 3 . Thus, the restriction 
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introduced by the kinematical boundary of the phase space for the two-body decay with the 
modifications of the non- logarithmic terms for Process (G2) given in Ref. 4) cannot be used 
in the region satisfying z 3 <^Z2 < ^z\ for the parton shower model. 

In order to cancel the negativity of the contribution from the modified three-body decay 
function [Eq. (5- 12a)] in the soft anti-quark region for Process (G2), we should construct 
another method. Since the negativity of the modified three-body decay function comes from 
the soft anti-quark region, we investigate the dominant term in the soft anti-quark region. 
In the soft anti-quark z% region, the logarithmic term Vl is given by Eq. (4-8), and the terms 
proportional to each of the logarithmic terms in Vl can also be approximated by 

L L 1) ^C F T R H{y 2 )-, (5-14) 
Lf~2C A T R —, (5-15) 

Z2V2 

I®~-C A T R —. (5-16) 

Z2V2 

From Eqs. (5-15) and (5-16), the logarithmic term proportional to the color factor C A T R is 

1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 



-0.2 
-0.4 

10" 4 10" 3 10" 2 10" 1 

Fig. 14. The Z2 dependences of Eq. (542a), taken from Ref. 4) and Eq. (5-12a) using Eq. (340) 
for Process (G2) with z 3 = 10~ 3 . Here, z 2 =10~ 4 -0.9. 

given by 



Eq.(5.12a) using (3.10) 
Eq.(5.12a),taken from Ref.4) 




R- 



Z2V2 



■ log Zi - CaTr — - — log 



Z2V2 



Z\Z 2 
Z 3 



= 2C A T R —\og( 



Z2V2 



ziz 3 y/2 

Z 2 ' 



(5-17) 
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This term Eq. (5-17) can be absorbed into the kinematical boundary of the phase space for 
two-body decay as follows: 



f d -hl2C A T R ^ + 2C A T R ^\ 0g (^) 112 

JMl t 13 Z 2 %)2 Z 2 %)2 V Z 2 J 

J 

J A 



,; " 3)1/2rftl3 2^T^. (5-18) 



'm 2 t 13 z 2 y 2 

After the absorption of the most singular term for the soft anti-quark region into the phase 
space of the two-body decay, the restriction on the phase space reduces to £13 < t{z\z 3 / Z2) 1 ^ 2 . 

In Fig. 15, since the contribution of the logarithmic term is small in the soft anti-quark 
region, the positivity of the three-body decay function for Process (G2) cannot be recovered 
by the subtraction of the effect of Eq. (5-17). We need other modifications. As explained 
below, by the use of Eq. (5-11), the logarithmic terms and the convolution terms of the non- 
logarithmic terms for the ladder diagrams and can be absorbed into the kinematical 
boundaries of the two-body decays. 

As explained above, we examine the contribution of the three-body decay function for 
the Process (G2) modified by Eqs. (5-11) and (5-18). Here, the ratio R qgq - is defined as 



log Vl + Lf logy, + (Jjg + L$)(CV) 

Rqgq = T} • O iy J 

VLL 

Using Eq. (5-19), the modified three-body decay function is given as 



(z 2 y2)V LL z 2 

where R qgq - is given by Eq. (5-19) and R\ is given by Eq. (4-1). In Fig. 15, the contribution of 
Eq. (5-20) is positive for small z 3 . Equation (5-20) represents the three-body decay function 
modified by the new restriction resulting from the kinematical boundary of the phase space 
for the two-body decay. The function in Eq. (5-20) can be applied for the region satisfying 
-22^-23- Thus, the new restriction due to the kinematical boundary of the phase space for 
the two-body decay can be applied to the region in which the contribution of Eq. (5 -12a) is 
negative. 

5.3. Modified non-logarithmic term for Process (Ql) 

As shown in Fig. 4, in case that the gluon momentum fraction is hard in the soft quark 
region for Process (Ql), the effect of the soft gluon is small. Also, the contribution of the 
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10" 1 



Eq.(5.1 2a) using Eq. (3.10) + 
Eq.(5.20) x 
Eq.(5.1 7) subtracted * *y 

X x 
X ^ 



"++++- 




10- 



10" : 



10" ; 



10"' 



Fig. 15. The z 2 dependence of Eq. (5- 12a) using Eq. (3-10), that of Eq. (5-20), and that of 
contribution of the three-body decay function with Eq. (5-17) subtracted at z 3 = 10~ 4 for 
Process (G2). 

non-logarithmic term Vn gives an important contribution in the soft quark region. We should 
evaluate the contribution of the NLL order for the parton shower model using the corrected 
method. 

Following the approach used in Ref. 4) and employing Eq. (5-11), we adopt a modi- 
fied non-logarithmic term and obtain a positive value in the soft quark region. The non- 
logarithmic terms of the ladder diagrams for Process (Ql) are given as follows: 



Cl 



r(l) 
N 



vl 



m 



(i) 

J N 5 



'(3) 
J N 



2C pC a 



-P(z 3 )k(^)- + z 3 



2/3 ' 2/3 



(*i - z 2 f 

yi 



(5-21) 

(5-22) 
(5-23) 



Hence, the convolution terms of the Altarelli-Parisi splitting functions n ) in Eqs. (5-21)- 
(5-23) can be absorbed into the kinematical boundaries of the phase spaces of the two-body 
decays. ^ The modified three-body decay function is defined as 



v nll -Y.l { ${cv) 



V qGG 



i=l 



V, 



LL 



(5-24a) 



where 

3 



5>#(CV) = -C F P( yi )P(^)--C F P(y 2 )P(^)--2C F C A P(z 3 )K(^)-. (5-24b) 

~i v yi ' yi v 2/2 y 2/2 v 2/3 y 2/3 
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0.4 - 




X 



- 



(A.16) in Ref.4) 



-0.2 



Eq.(5.24a) 



X 



X 



-0.4 - 



-0.6 - 



X 



-0.8 



Fig. 16. The z 2 dependences of Eqs. (5-24a) and (A16a) taken from Ref. 4) with z% = 10 1 for 
Process (Ql). 

As shown in Fig. 16, the contribution of the modified three-body decay function [Eq. (5-24a)] 
is positive in the soft quark region. Although the correction of NLL order is larger than that 
reported in Ref. 4), the dominant contribution corresponding to the effect of the soft quark 
can be subtracted exactly. 



We have calculated three-body decay functions using the time-like jet calculus. Although 
most of the results agree with previous results obtained by two groups, some of the calculated 
results in this paper are different. The results we found for Processes (Gl) and (Q2) agree 
with those reported in Refs. 4) and 5). Also, the result for Process (Ql) obtained in this 
paper is consistent with that of Ref. 4). However, the results we found for 1^ and 1^ of 
Process (G2) are not equivalent to those given in Refs. 4) and 5). 

We also studied the properties of the three-body decay functions. First, we examined 
the gluon momentum fraction z 2 dependences of the ratios R\ and R 2 . We found that the 
contribution of the logarithmic term V L is dominant in most regions, with the exceptions of 
the soft quark region for Process (Ql) and the soft anti-quark region for Process (G2). 

Using numerical results, we examined the structure of the strong cancellation for the 
non-logarithmic term V N . Similarly to the situation considered in Ref. 7), we found that 
the contributions of the non- logarithmic terms for the interference diagrams are cancelled 
by those for the ladder diagrams and the other interference diagrams. 

We also studied the properties of the three-body decay functions in the soft parton 



6. Summary 
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regions. In the soft quark region for Process (Ql) and the soft anti-quark region for Process 
(G2), the contribution of the logarithmic term Vj, is small, because the contributions of the 
ladder diagrams and the interference diagrams cancel. 

Because some of our results for the three-body decay functions differ from the results 
given in Refs. 4) and 5), the contribution of the modified three-body decay function for 
Process (G2) derived in Ref. 4) was examined. We also corrected some of the results for the 
convolution terms of the non-logarithmic terms absorbed into the kinematical boundaries 
of the two-body decays given in Ref. 4). It was found that using the new results, the 
contribution of the modified three-body decay function used in Ref. 4) is negative in some 
regions. Therefore, we suggested a new restriction introduced by the kinematical boundary 
of the phase space for the two-body decay with modifications of the non-logarithmic terms. 
The result of this change is that the contribution of the three-body decay function for Process 
(G2) is positive in the region satisfying z^<^Z2<^z\. 

The contribution of NLL order in the soft quark region for Process (Ql) was examined. 
The non-logarithmic term Vn corresponding to the dominant contribution in the soft quark 
region was subtracted as the effect of the soft gluon in NLLJET. Considering this, we sug- 
gested that some of the non- logarithmic terms for the ladder diagrams should be absorbed 
into the kinematical boundaries of the two-body decays. 

To summarize, we derived new three-body decay functions for time-like branching and 
carried out a numerical treatment of the three-body decay functions in the soft parton regions 
for the parton shower model. In a future paper, we will present the NLL parton shower model 
improved by the new results and the methods presented in this paper. 
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Appendix A 

Kinematics 

In this appendix, we present the details of the calculations for the decay functions of 
Types [LI] and [13]. The calculations for the other types can be obtained by exchanging the 
momenta U in Types [LI] and [13]. 
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A.l. Calculation of Type [LI] (Ladder diagram) 

The combination X = Z\Z 2 is used to determine the decay function for Type [LI]. The 
constraint in the 5(1%) function is written 

= ll = (q ~ k ~ hf = z 3 t + (-^-)l 2 1T - V -V\ T . (A-l) 

v Z\V\> z 2 

— * — * 

Here, 1 2 t has been replaced by introducing the new vectors p 2 r and 1\t according to 



l 2T = -—l 1T +p 2T , (A-2) 

yi 



with yi = 1 - Zi (i= 1,2,3). 

From Eq. (A-l), p\ T is given by 



Z2Z3 ( . I 2 



I \ T is written as the invariant 



l\ T = ztyit - zitw (A-4) 
By use of Eq. (A-4), p\ T is given by 

PIt = (A-5) 

It is easy to calculate the term P^M in Eq. (2-9) using the constraint represented by the 
delta function. The general form of the term P^M is given by 

P&M — + -^12 + -^3^23 + -^4^12^23 + -^5^12 + -^6^23 (A-6) 

^23 

where the Fj are functions of the Z{. The variable t\ 2 depends on the azimuthal angle <fi. The 
variable t i2 is replaced by Eq. (A-2): 



/, t \9 ( Z 2 Z\Z 2 2,Z 2 \->n Z\ _,n 2 , , , r» , 

t\i = (h + h) = 1 2~ H )ti T H P2T p 2 T Kit cos0. 

j/i yi ' z 2 Vl 

(A-7) 
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Here, the term with dependence on the azimuthal angle in Eq. (A- 7) can be isolated. It 
follows from Eqs. (A-4) and (A-5) that tu and t\ 2 in the numerator satisfy 



z 2 Z\Z 3 — Z 2 2 . 

*i2 = — t H 72 ~ —\P2t\\Iit\ cos0, (A-8) 



t 2 

L 12 



yi v\ yi 

z 2 Z1Z3 — Z 2 ]2 i A Z\Z 2 Z 3 
1 H 9 ^23 

y\ yi 



+ -t 23 )t 23 cos 2 0. (A-9) 

y\ v 7 



After replacement of the invariant in Eq. (A-6), we obtain 



p(i)M = Cit 2 + G 2 t 2 23 + G 3 tt 23 + {G A t 2 23 + gsM cos 2 + (G 6 t 2 3 + G 7 t) cos 

(A-10) 



t 23 



where the Gj are functions of the Zj. Although G 6 and G7, which are proportional to cos0, 
are also functions of f>2T and /it, they drop out through the azimuthal integration. Since 
the mass singularity is of logarithmic order, there is no singularity term of higher than 
logarithmic order. Thus, we obtain the condition 

Gi = 0. (A-ll) 

Next, using Eqs. (A-4) and (A-5), the phase space volume element is given by 

dt = -zxdt^dpl^dill). (A-12) 
2 

The integration of the 8 function is trivial. The kinematical boundary of the integration 
over the variable t 23 is determined from the conditions p \ T > and I \ T > 0. We obtain 

< t 23 < Vl t. (A-13) 
Therefore, the formula for the phase space volume element in the jet calculus is given by 



dr 



Z\Z 2 ix 2 pi* d(f> 



yi Jul 2-k 

We integrate the term P®M over the azimuthal angle to obtain 



f d Ap&M= H2t ^ + 2 H ' tt2 \ (A-15) 
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where the Hi are functions of the Z{. After the integration over the azimuthal angle, the 
terms proportional to cos 2 in Eq. (A-10) are absorbed into the other terms. In order to 
obtain the decay function J in Eq. (2-19), integration over the variable £23 is performed. We 
then obtain 




(A-16) 



The first term and the other terms in Eq. (A-16) correspond to the contributions of LL order 
and NLL order, respectively. The contribution of NLL order is written 

-^-H,\o gyi + \h 2 = logj/x + L%\ (A-17) 

with 

lP = ±H 3 ,L$ = \h 2 . (A-18) 

Since the two-body decay function includes 0(M%/t), we do not need this term in the three- 
body decay function. 



A. 2. Calculation of Type [13] {Interference diagram) 

The decay functions for the interference diagrams are free from mass singularities. Here 
we construct the condition for these decay function to possess no mass singularity. 
We start by constructing the phase space volume element as follows: 

dP = -d(l\ T )d(ll T )d<p5(ll). (A- 19) 

— * — * 

The new vectors l 1T and 1 2 t are defined as 

hir = kr, hi T H h 2 T = kr- (A-20) 

V\ ZiVi 

We find the phase space volume element from Eq. (A-20): 
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dP = l(^) 2 d(h 2 1T )d<t>d(hi T )6(ll). (A-21) 



The constraint represented by the 5 function is expressed as 



>2 

ziyi " z 2 zfyi 



= l 2 3 = z 3 t - ^-h\ T - ^r-h\ T . (A-22) 



Next, we calculate the term P^M in Eq. (2-9). After eliminating t\ 2 using Eq. (2-8), 
the general formula for the term P^M becomes 

p& M = Flt2 + -^2^13^23 + -^3^13 + -^4^23 + -^13 + -^23 (A-23) 

^23^13 

We absorb the last two terms into the other terms in Eq. (A-23). The last two terms are 
given by 



K 1 =fdT^- = fdP^ (A-24) 

J t ; 13' : 23 J *23 



and 



" ^23 . fj^3 



k 2 = 1 dr-^- = i dr-^. (A-25) 

J ^13^23 J ^13 

After z 2 in Eq. (A-8) is replaced by z$, the invariant £i 3 is given by 

ti 3 = — t + 9 Z3 t 23 - —\p3t\\Iit\ cos (p. (A-26) 
V\ Vi Vi 

Similarly, z\ in Eq. (A-26) can be replaced by z 2 . The invariant £23 is given by 

, Z3 z\z 2 — Z3 2 . „m 

^23 = — t H 2 *13 \P3T «2T COS 0. (A-27) 

1/2 2/2 2/2 

Using Eqs. (A-26) and (A-27), we obtain the formula for the term P®M as 



p«M = J dP 



ft £1*13*23 + G 2t 2 + G 3 ttl3 + ^4*23* + (G 5 t23 + G 6 t 13 ) COS 



*13*23 

(A-28) 



where the Gi are functions of the Z{. G5 and G% are functions of p 3 T and lyr, 1 2 t- The terms 
proportional to cos can be dropped due to the azimuthal integration. 
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The interference term must be free from mass singularities, which may occur at t\3 = 
or ^23 = 0. The term G\t 1^23 in Eq. (A-28) has no mass singularity. The terms proportional 
to G2, G3 and G4 must also be free of mass singularities. Here we present a non-trivial 
combination of the functions G; which results in no mass singularity. The invariants t±3 and 
t 2 3 are written 



1-23 



2„ 2Ti 



Z\Z2 



(A-29) 



Z 3 (7 



'-13 



(K 1T + h 2T y 



(A-30) 



Using Eqs. (A-4), (A-29) and (A-30), the formula for the azimuthal integration satisfying 
the condition that there be no mass singularity is given by 



2yizfz 2 f 2w h\ T + h 1T h 2T 



rZTT 

/ d<j> 
Jo 



h 2 2T (hiT + h 2T f 



^V\ Z \ Z 2 al t2 



■25 



6(h 2T — h 1T ) 



/ dv—r-yt *ia ^23, 



(A-31) 



r-2vr 

'0 fl3 r 23 Z3 Z3 

where 9 is the step function. The relation between Eq. (A-28) and Eq. (A-31) is given by 



f~i 2/l f~i f~i 1/2 f~i 

\J3 — Or 2, <-T4 — U-2 

Z3 Z3 

Substituting this combination into Eq. (A-28), we have 



(A-32) 



[dr\G 1 + G2-^-(t- ^3-^23) 

J ^13^23 ^3 Z3 



Then, the latter integrations can be performed: 



(A-33) 



K = 
K 3 = 



,f^l3^23 
^13^23 

[df—(t 



yi, V2. \ 

— tl3 ^23 ) • 

Z3 Z3 > 



(A-34) 
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The P®M term in Kq is independent of the azimuthal angle, while the P^M term in A3 
depends on the azimuthal angle. To begin, we integrate the phase space in Kq. Substituting 
Eq. (A-29) into Eq. (A-21) and integrating over hf T , we have 



df=^ (A-35) 
Hi Jo 2n 

The region of integration over the invariant t 2 3 is the same as that in Eq. (A- 13). We obtain 



K = ir 2 Zl z 2 t. (A-36) 

Next, we compute A3. The result obtained from integrating A3 over the azimuthal angle 
corresponds to the second term in Eq. (A-31). Due to the constraint introduced by the S 
function and the step function [Eqs. (A-22) and (A-31)], the kinematical boundary of A 3 is 
given by 



h 2 1T >0^t 23 < Vl t, (A-37) 

V2 



h 2 1T <h 2 2T ^—t<t 23 . (A-38) 



Therefore, the result for K 3 is 



A 3 = -**™t r* — = -^^tlog^. (A-39) 

Z "i J%t t 2 3 Z3 Z3 



Thus, the decay function of Type [13] is given by 



A-K 2 ZiZ 2 t 

ic.-J-G, log m, (A40) 

4 2^3 z 3 



with 
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